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EOMET  '77 

Atmospheric  Transmission  Values  Derived  from 
Photographs  of  the  Sea  Horizon 

1.  INTRODUCTION 

1. 1 There  is  limited  experimental  evidence  that  the  transmission 

of  the  atmosphere  at  sea  during  the  daytime  can  be  deduced  from  photo- 
graphs of  the  horizon.  ' 2 The  experiments  which  have  been  reported 
result  in  reasonable  values  of  transmission, but  no  separate  measurement 
of  transmission  or  transmission-related  values  have  been  made  con- 
current with  the  horizon  photographs  to  demonstrate  that  the  technique  is 
valid. 


1.2  The  1977  EOMET  Cruise  presented  an  opportunity  to  compare 
the  results  of  horizon  photographs  with  concurrently  measured  or  re- 
corded meteorological  data  such  as  visibility,  aerosol  concentration  and 
size  distribution,  nephelometer  readings  etc.  To  test  the  techniques  and 
implement  such  comparisons,  groups  of  horizon  photographs  were  made 
on  twenty  occasions  between  May  17  and  June  6,  1977. 

1.3  One  photograph  from  each  of  the  twenty  groups  has  been  se- 
lected by  visual  inspection  to  be  most  promising  for  data  reduction.  In 
addition,  a properly  exposed  but  challenging  photograph  has  been  selected 

from  ten  of  the  groups.  Challenge  includes  cloud  shadows,  wind  streaks, 
and  white  caps. 


1.4  The  horizon  photographs  were  made  with  color  filmland  data 

reduction  on  a particular  film  results  in  three  values  of  atmospheric 


’» 


scattering  coefficient,  one  in  the  blue  centering  on  X = 0.  44p,  one 
in  the  green  with  X = 0.  55p  and  one  in  the  red  at  X = 0.  64p.  The 
coefficients  are  appropriate  for  the  sum  of  molecular  and  aerosol  scat- 
tering. Since  molecular  coefficients  are  known  the  aerosol  scattering 
coefficients  can  be  obtained  if  the  method  is  valid. 

1.5  The  end  product  of  this  report  is  a tabulation  of  aerosol 

scattering  coefficients  at  three  wavelengths  obtained  from  color  photo- 
graphs of  the  horizon  taken  on  twenty  occasions  during  the  1977  cruise 
of  the  USNS  Hayes.  In  addition,  challenging  photographs  from  ten 
occasions  have  been  analyzed  and  the  scattering  coefficients  obtained  from 
them  are  compared  with  coefficients  obtained  from  the  corresponding 
most  promising  photographs.  Pertinent  supporting  data  are  also  presented. 
This  includes  date,  time  of  day,  latitude,  longitude,  wind  direction, 
ship  heading,  and  cloud  type  and  cover. 
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2. 


THEORY 


2. 1 Equations  for  the  apparent  radiance  of  the  sea  at  and  near 
the  horizon  are  developed  here.  The  equations  show  the  relationship 
between  apparent  radiance  and  atmospheric  scattering  coefficient.  In 
particular,  if  photographs  of  the  sea  horizon  are  taken  from  a known 
altitude  above  sea  level  then  relative  radiance  values  extracted  from 
the  photographs  may  be  used  to  calculate  scattering  coefficients. 

2.2  Figure  1 illustrates  the  observing  conditions.  In  the  figure, 

P is  the  radius  of  the  earth,  R is  the  distance  to  the  point  observed  on  the 
sea  surface,  h is  the  height  of  the  observer,  9 is  the  depression  from  the 
horizontal  of  the  line  to  the  horizon  and  9 + <t>  is  the  depression  from  the 
horizontal  of  the  point  observed.  In  what  follows  equations  are  developed 
for  the  apparent  radiance  of  the  path  from  the  observer  to  any  point  on 
the  sea  surface.  These  equations  are  so  written  that  if  the  radiance  is 
measured  for  various  values  of  0 then  a value  for  the  scattering  coefficient 
can  be  obtained.  Following  this,  equations  relating  R to  0 are  developed. 
Finally  the  relationships  between  positions  on  the  picture  of  the  horizon 
and  values  of  0 are  discussed  as  well  as  the  relationship  of  exposure  of 
the  film  at  those  positions  to  apparent  radiance  of  the  sea  surface. 


2.3  In  equation  1 which  follows,  first  order  scattering  theory  is 

used  to  describe  the  apparent  radiance  N (0)  of  the  sea  surface  at  angle 


0 below  the  line  to  the  horizon. 


N(0)  = 


(x(r)dr 


rR 

o<r)dr 


E(r)  <3>(r)o(r)  e 


dr+E(r)g(r)e 


I 


where  N (0)  = the  apparent  radiance  of  the  sea  surface  at  angle 

0 below  the  line  to  the  horizon. 
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Geometry  used: 

h ■ height  observer's  eye  in  meters 
0 ~ dip  of  the  horizon  in  minutes 
9 ++  - dip  of  point  in  question  in  minutes 
R = distance  to  point  in  kilometers 
From  Bowditch: 

982  (9  ++)*.  15<35  h 


R 


2.232  (6  + * ) 


E (r)  = total  irradiance  at  position  r along  path 
ct  (r)  = the  effective  phase  function  for  scattering 
from  position  r toward  observing  point. 
a (r)  = scattering  coefficient  at  r. 
g (R)  = a coefficient  such  that  g (R)  E (R)  is  the 

inherent  radiance  of  the  sea  surface  in  the 
direction  of  the  observing  point. 


2-4  Equation  1 can  be  solved  if  three  assumptions  are  made. 

The  first  is  that  the  irradiance  along  the  path  is  constant,  the  second  is 
that  the  phase  function  along  the  path  is  constant,  and  the  third  is  that 
the  inherent  radiance  of  the  sea  surface  is  the  same  for  all  positions  ob- 
served. With  these  assumptions  equation  1 may  be  written 

r)dr  - fE  CT(r)  dr 

Jo  Jo 

N(*)  = E<J>[1  -e  ] + E g e <2) 


The  derivative  of  this  radiance  as  a function  of  the  range  to  the  sea 
surface,  R,  is 

_J cr(r)dr 

^ N (*)  - E (<J>  - P ) g (R)  e 0 


When  the  data  has  been  taken  in  the  presence  of  fog  there  is  some  range 

^ 3$C 

R (corresponding  to  some  value  of  * = <p  ) over  which  the  transmission 
is  essentially  zero.  This  implies 


In  clear  weather  the  horizon  can  be  seen  so  there  is  no  path  to  the  sea 
surface  satisfying  equation  4.  However,  the  path  just  above  the  horizon 
is  assumed  to  involve  sufficient  atmosphere  that  equation  4 holds.  For 
these  two  cases,  represented  by  0*  which  is  positive  in  fog  and  negative 
in  clear  weather,  substitution  of  equation  4 in  equation  2 leads  to 


N (0  ) = E0 


The  difference  between  equation  5 and  equation  2 is  given 


in  equation  6 


N (0  ) - N (0)  = E ( <J>  - p ) e 


o<r)dr 

~ Jo 


Finally,  equation  3 divided  by  equation  6 gives  the  value  of  the  scattering 
coefficient  at  distance  R.  This  is  shown  in  equation  7 


tr(R)  = 


N (0*)  - N (0) 


N (0) 


This  may  also  be  written 


cr(R)  = - 


N(0r)  - N (0)  {N(0*)-N(0)J 

= 'di~  ln  {N  <**>  “ N (♦>} 

= - ^ [In  |n  ( 0*)  - N (0)  | - ln  N ( 0*)  ] 


= d ^ JNJf  ) - N(0)  \ 
dR  \ N (0*)  / 


> / 


This  is  the  key  equation  for  the  work  discussed  here.  Further  dis- 
cussion of  equation  8 and  its  application  to  the  data  reduction  will  be 
made  easier  by  an  understanding  of  the  relationship  of  R to  0.  Equations 
showing  that  relationship  are  developed  next. 

2.  6 Two  methods  are  given  here  for  calculating  the  distance  to  a 

point  short  of  the  horizon  when  the  height  of  the  observed  is  known  and 
the  angle  of  depression  of  the  point  relative  to  the  horizon  is  measured. 

The  first  method  has  its  origin  in  geometrical  calculations  with  no  cor- 
rections for  atmospheric  refraction  and  the  second  method  used  empirical 
relationships  which  are  sensitive  to  refractive  effects  and  are  taken  from 
Bowdich^. 

2.7  In  Figure  1 an  observer  at  O views  at  distance  R a point  p short 

of  the  horizon  b.  The  dip  angle  of  the  horizon  is  0 and  the  dip  angle  of  the 
point  is  0 + $.  The  altitude  of  the  observer  is  h and  the  radius  of  the  earth 
is  p . Using  the  cosine  law  one  may  write 

+ ( p + h)2  + R2  -2(  P + h ) r cos(  W 2 - ( 0 + $)  ) 

2 2 
or  R-2(p+h)rsin(0  + $)+2ph+h=O  (9) 

Solving  for  R gives 

R = ( p + h ) sin  ( 0 + <t>  ) + y(p+h)2sin2(0+*)-2ph-h2  (10) 

I 

The  solution  with  the  positive  radical  is  appropriate  to  point  p in  Figure 
1.  The  solution  with  the  negative  radical  is  appropriate  for  this  develop- 
ment as  it  relates  to  the  point  p which  is  short  of  the  horizon.  When  the 

f 

radical  is  equal  to  zero  the  points  p and  p are  coincident  and  the  line  of 


10 


sight  is  directed  to  the  horizon.  In  this  case 


( P + h)2  sin2  (0  + 0)-2ph-h2  = 


2 2 ph  - h 

°r  Sine=(.  + h)2 


-1  2ph  ~ h 
and  0 ~ sin  \ 2 

I ( P + h) 


This  is  the  dip  of  horizon  assuming  no  refraction.  Now  equation  10 
may  be  written 

./  sin20 

R = ( P + h)  sin  ( 0 + 0)  \1  - .2, 

1 sin  ( 0 + 6 ) 


Equations  11  and  12  make  up  the  geometrical  solution  for  the  relationship 
between  R and  0. 

2.  8 Bowditch  gives  an  expression  for  the  dip  of  a point  short  of  the 

horizon  which  is  rewritten  here  in  metric  units  as 

(0  + 0)  = 0.224R  + 3.440  h/R  (13; 

where 

( 0 + 0 ) is  the  dip  of  the  point  in  minutes 
h is  the  height  of  the  observers  eye  in  minutes  . 

and  R is  the  distance  to  the  point  in  kilometers. 

Solving  equation  13  for  R,  the  distance  to  the  point  short  of  the  horizon 
gives 

R = 2.232  (0  +0)  - 1/4. 982  (0  +0)2  - 15.35h  (14) 


' tpy, •r' */  • 


iisO 


For  the  case  0 = 0 (the  line  to  the  horizon) 


r 4 . 9826  - 15.  35  h = 0 


6 - 1.76  ^h~ 


2.  9 For  the  1977  EOMET  work  the  value  of  h,  the  height  of  the 

camera  lens  above  the  sea  surface,  was  4.  60  meters.  Usingthe  geo- 
metrical equations  gives  0,  the  depression  of  the  horizon,  as  4.  13 
minutes  and  the  distance  to  the  horizon  as  7,656  km.  Using  Bowditch's 
equations  gives  0 as  3.  775  min.  and  the  distance  to  the  horizon  as  7.  810  km. 
In  Figure  2 the  range  to  the  horizon  calculated  using  Bowditch's  equation, 
RBowd’  1S  shown  as  a function  of  0 . The  figure  also  shows  the  ratio  8 where 


Bowd 

R_ 


Geom 


Bowd 

In  the  data  reduction  done  here  Bowditch's  equations  are  used. 

2. 10  Equation  8 relates  cr(R),  the  scattering  coefficient  at  range 

R,  to  the  observable  radiances  N (0*)  and  N (0o)  where  N ( 0*)  is  the 
radiance  just  above  the  horizon  or,  in  fog,  at  a dip  angle  0*  small  enough 
so  that  the  specular  transmission  to  the  sea  surface  approaches  zero. 
Equation  8 is  repeated  here 


o(R)  = 


. ( N(0*)  - N (0)  ) 
( N (0*)  ) 


Equation  14  is  the  expression  from  Bowditch  forR,  km,  the  range  to  a 
point  short  of  the  horizon,  when  the  point  has  a dip  ( 6 + 0 ),  min.  the 
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Figure  2.  Comparison  of  Geometric  Solution  to  Bowditch 


horizon  has  a dip  6,  min  and  the  observer  is  a height  h,  m,  above  the 
sea  surface.  As  shown  in  Figure  2.  when  h = 4.  6 m the  relationship 

°f  R to  ♦ *n  *he  range  1 < $ < 10  min  can  be  given  by  the  approximate 
expressions 


R ~4.1  +3  log  $ 


(16a) 


and 


so 


dR  = 
d$ 

J_  dR 
R d$ 


1.3 

$ 

-1.3 

$(4.1  - 3 log  $) 


(16b) 

(16c) 


2. 11  In  the  EOMET  experiments  h was  4.  6 m and  the  lens  of  the 

camera  used  had  a focal  length  of  50mm.  The  resolution  of  the  film  used 
was.  at  best,  50  lines  per  mm.  In  a rough  way  the  angular  resolution  may 
be  taken  as  A a - tan  — = 1.  38  min.  The  probable  error  in 
assigning  location  of  the  horizon  is  of  the  class  of  ± one  resolution  ele- 
ment, or  ± 1.38  min.  Table  1 presents  calculations  using  equation  16c 
of  A R the  uncertainty  in  range  R to  points  short  of  the  horizon  due  to 
+ A a , the  uncertainties  in  the  location  of  the  horizon.  One  film  resol- 
ution element  is  taken  as  the  value  of  Aa . Values  of  A R are  calculated 
for  true  distances  on  the  film  below  the  horizon  from  one  to  ten  resolution 
elements  in  length. 


2.12  In  the  EOMET  experiments  there  was  no  effort  made  to  hold 

the  camera  exactly  horizontal  or  to  point  it  exactly  at  the  horizon.  Location 
of  the  horizon  image  on  the  films  introduces  uncertainties  in  angular  posi- 
tion of  points  short  of  the  horizon  of  the  order  of  + 1.  375  min.  This  un- 
certainty makes  the  determination  of  scattering  coefficient  as  a function  of 


'itgrer 


TABLE  1.  ERROR  ANALYSIS 


Number  of  film 
resolution  elements 
from  horizon 

<t>,  min 

R,km 

A 

ar 

R 

0 

0 

7.810 

1 

1.  375 

3.  685 

+ 

1.3 

+ 0.  35 

2 

2.750 

2.782 

± 

0.  65 

t 0.23 

3 

4.  125 

2.254 

+ 

0.43 

± 0.  19 

4 

5.  500 

1.879 

. t 

0.33 

t 0.17 

5 

6.875 

1.  588 

± 

0.  26 

± 0.16 

6 

8.250 

1.  351 

+ 

0.22 

+ 0.  16 

7 

9.  625 

1.150 

+ 

0.  19 

+ 0.16 

In  this  table  AR  is  the  uncertainty  in  range  R to  point 
short  of  horizon.  True  point  position,  $,  in  minutes,  is  asso- 
ciated with  the  number  of  resolution  elements  on  film  below 
the  horizon  image.  The  uncertainty  in  locating  the  horizon  is 
± Aq.  where  A a is  taken  to  be  equal  to  one  resolution  element 
which  for  the  lens  and  film  used  is  1.  375  min. 
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range,  <t(R),  by  use  of  equation  8 essentially  impossible  In  a fall-back' 
position  it  is  assumed  that  the  scattering  coefficient  a(R)  = a;  that  is, 
that  the  scattering  coefficient  is  independent  of  range  then,  from  equation 
8 


Ih  ■«>  = 

dR 


( N (0*)  - N (0)  ) 
\ N (**)  ) 


and 


d i 

— In 
dR 


( N (»*)  - N (»)  ) 
( N (**)  ) 


= const  * a 


or  a plot  of 
- In 


fN  N(»)  I VSR 

( N («  ) ) 


(17) 


(18) 


must  be  a straight  line  with  slope  equal  to  a 
2.13  Figure  3 presents  plots  of  the  function 

f(*)  = In. 


(N(»*)  - N(») ) 
I N ($*)  j 


vs 


R ( 0 + A*) 


(19) 


where  <f>  is  the  correct  value  of  the  angle  between  the  line  to  the  point  in 
question  and  the  line  to  the  horizon  and  is  the  error  in  that  value. 
Curve  A has  no  error  and  is  a straight  line  as  required  by  equation  18. 
Curves  B have  errors  of  ± 0.  34  min  (corresponding  to  0.  25  pixels  on  film) 


16 


and  curves  C have  errors  of  ± 0.  69  min.  In  this  illustration  the  range  to 
the  horizon  is  7.  97  km  and  the  correct  value  of  R for  point  1 is  4.  62  km. 
There  is  nothing  unique  about  this.  The  positive  or  negative  curvature  of 
the  plot  which  results  from  incorrect  guess  of  horizon  image  position  is 
characteristic  of  any  case  of  R < 7.  97  km.  This  feature  will  be  dis- 
cussed again  when  the  techniques  of  data  reduction  are  described. 
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3. 


DATA  REDUCTION 


3. 1 Ektachrome-X  film  was  used  for  the  EOMET  horizon  photo- 

graphs. All  of  the  film  had  the  same  emulsion  number  and  it  was  all 
processed  by  Eastman  Kodak  in  a single  batch  along  with  appropriate 
sensitometric  exposures.  Ektachrome-X  is  a positive  color  film  and  by 
using  appropriate  optical  filters  in  data  reduction  the  film  response  in 
three  wavelength  bands  can  be  determined.  These  bands  center  on  wave- 
lengths 0.44ju,  0.  55^,  and  0.  64p.  A Photometries  Data  Systems  Micro- 
densitometer was  used  in  data  reduction  and  the  first  step  of  the  process 
was  to  use  the  sensitometric  exposures  to  establish  at  each  of  the  three 
wavelengths  the  relationship  between  film  density  and  relative  exposure. 

3.  2 Figure  4 is  a sketch  showing  the  path  on  the  processed  horizon 

photographs  scanned  by  the  microdensitometer.  The  densitometer  slit, 
which  wat»  1300^  long  and  17p  wide  was  set  parallel  to  the  horizon  and 
moved  from  a position  one  millimeter  above  the  best  guess  of  horizon 
image  location  to  a position  5 mm  below  that  location.  The  density  of 
the  film  was  read  automatically  every  20n  along  the  track  and  these 
readings  were  recorded  on  tape.  This  process  was  carried  out  for  each 
of  the  three  colors,  blue,  green,  and  red  which  were  selected  by  incorpo- 
rating in  the  optical  train  of  the  densitometer  Wratten  filters  No.  94, 

No.  93,  and  No.  92  respectively. 

3.3  A simple  program  was  set  up  for  a CDC  6600  computer 

which  used  the  sensitometer  data  to  convert  density  to  relative  exposure 
for  each  of  the  three  colors.  The  product  of  this  program  is  a properly 
annotated  print  out  of  the  step  number  and  relative  exposure  associated 
with  the  microdensitometer  track  across  the  horizon  image.  Figure  5 
is  such  a print  out  and  is  for  frame  13  of  the  film  exposed  on  24  May  i977 
at  1355.  This  microdensitometer  run  was  made  using  a Wratten  Filter 
No.  94,  which  is  blue.  The  best  guess  of  horizon  position  was  assigned 

. \ 
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Figure  5.-  Values  of  relative  exposure  from  horizon  *race; 
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track  distance  scale  value  of  10.  00  mm.  The  trace  was  started  above 
the  horizon  image  at  9.  00  mm  and  run  to  15.  00  mm. 

3.4  The  relative  exposure  values  for  the  15  points  between 

9.  60  mm  and  10.  84  mm  were  averaged  (10.  89)  and  used  as  N (0*)  the 
relative  radiance  of  the  «sky  above  the  horizon.  The  exposures  in  the 
sequence  starting  at  9.  98  were  scanned  visually  and  sequential  differences 
noted.  The  exposure  (8.  68  at  10.  04)  following  the  greatest  sequential 
difference  (10.  37-8.  68)  was  selected  as  the  first  of  a series  of  nine 
values  of  N (0)  to  be  used  in  calculating  the  scattering  coefficient. 

5 A program  for  an  HP-97  was  used  to  calculate 


f(x  - x ) 
o 


-In 


(n  (0*)  - N (0)  ) 
( N (0*)  ( 


(20) 


for  each  of  the  nine  values  of  N (0)  and  also  the  ranges,  R,  fdr  each  of 
the  nine  values  of  (x-x^)  in  the  microdensitometer  track  assuming  some 
position  xq  for  the  horizon.  The  HP-97  program  calculated  the  linear 
regression  corresponding  to 


f(x-XQ)  = A + BR  (x-Xq)  (21) 

for  the  first  five  values  of  (x-x  ).  yielding  B and  then  for  the  nine 

o 5 

values  of  (x-xq),  yielding  Bg.  The  ratio  B5/Bg  was  calculated  and  if 

this  was  within  one  percent  of  unity  the  selected  position  for  the  horizon 

was  accepted  otherwise  a new  position  for  the  horizon  (x  ) was  chosen 

o 

Usually  no  more  than  four  such  iterations  of  this  sort  were  necessary.  The 
final  value  of  Bg,  the  slope  of  the  regression,  was  taken  as  a,  the  scattering 
coefficient.  The  program  then  could  print  out  the  nine  pairs  of  values 
N (0)  and  R appropriate  to  the  nine  values  of  (x-xq).  The  following  values 
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■ iF# 


Jr 

- 


Tyr*. 


from  Figure  4 were  used  in  this  process: 


* 


N 

(♦  ) 

= 

10.89 

N 

(*>.  1 

= 

8.  68 

at 

10.04 

2 

= 

7.  53 

10.06 

3 

= 

7.  15 

10.08 

4 

* 

6.  78 

10.  10 

5 

* 

6.  78 

10.  12 

6 

= 

6.  61 

10.  14 

7 

= 

6.  61 

10.  16 

8 

= 

6.  35 

10.  18 

9 

= 

6.  19 

10.20 

In  the  iteration  to  find  the  horizon  position  the  following 

values  of  x were  used  resulting  in  the  associated  values  of  B /B 
° 5 9" 

x 

o 

10.03 

10.04 
10.  375 
10.039 
10.0395 


VB9 

1.  18 
0.975 
1.06 
1.02 
1.008 


lue  x = 
o 

10.0395  values  of  f(x- 

R,km 

f(x-x 

o 

7.  18 

1.  59 

3.  62 

1.18 

2.  65 

1.07 

2.  12 

0.974 

1.77 

0.974 

l.  53 

0.  934 

1.  35 

0.934 

1.20 

0.850 

1. 09 

0.  840 

The  regression  is  given  by 

f(x-x  ) = 0.744  + 0.119  R (x-x  ) (22) 

o o 

hence  a = 0.119  ~ 0.12  km  . 

3.6  In  this  work  30  films  have  been  analyzed  in  three  colors 

leading  to  90  possible  values  of  a.  Of  these  53  have  been  done  as 
described  above,  12  have  been  done  using  nine  averages  of  sequential 
groups  of  five  (45  values  of  N ($)  in  total),  3 have  been  done  with 
sequential  groups  of  fifteen  and  one  with  sequential  groups  of  3.  Of 
the  remaining  21,  18  have  been  done  by  hand  and  involved  less  than 
five  groups  of  from  1 to  15  values  of  N ($)  per  group  and  three  defied 
analysis.  Values  of  a so  obtained  are  presented  in  Table  2 . 
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TABLE  2.  Values  of  Scattering  Coefficient,  km 
Date  Time  Demonstration  Challenge 


0.  44m 

0.  55m 

0.  64m 

0.44^ 

0.  55m 

0.  64m 

May  17 

1320 

2.42 

1.95 

1.63  (5,9) 

18 

0945 

28 

8.  5 

3.5  (15,6) 

3.48 

3.  17 

2.79  (15,9) 

18 

1700 

0.  12 

0.11 

0.10  (1,9) 

0.  20 

0.  15 

0.08  (1,9) 

21 

1545 

0.  21 

0.  17 

0.10  (1,9) 

0.  36 

0.24 

0.28  (1,5) 

22 

1230 

0.  19 

0.22 

0.20  (1,9) 

0.  62 

0.  52 

0.51  (1,5) 

22 

1830 

0.08 

0.06 

0.08  (1,9) 

0.  23 

0.  29 

0.29  (1,9) 

23 

1100 

0.  16 

0.  13 

0.12  (1,9) 

24 

1030 

0.25 

0.  18 

0.19  (1,9) 

24 

1355 

0.  12 

0.08 

0.06  (1,9) 

1. 45 

1. 25 

1.29  (5,9) 

25 

1055 

0.  17 

0.  19 

0 17  (1,9) 

25 

1610 

0.  13 

0.085 

0.086(1,9) 

26 

0925 

0.27 

0.33 

0.21  (1,9) 

0.  51 

0,42 

0.45  (1,5) 

26 

1530 

0.43 

0.  31 

0.25  (1,9) 

0.41 

0.  36 

0.47  (1.4) 

27 

0950 

4.09 

2.97 

1.99  (5,9) 

31 

1050 

0.37 

0.29 

0.19  (1,9) 

Failure 

June  1 

1240 

0.  16 

0.14 

0.12  (1,9) 

3 

1230 

0.  16 

0.  15 

0.19  (1,9) 

4 

1505 

0.  56 

0.45 

0.35  (1,9) 

5 

1440 

0.  67 

0.  53 

0.38  (1,9) 

1.  61 

0.  92 

1.08  (5,4) 

6 

1330 

0.92 

0.91 

0.85  (5,9) 

DEMONSTRATION  VALUES  ARE  FROM  FRAME  SELECTED  VISUALLY 
AS  EASIEST  TO  REDUCE  . 

CHALLENGE  VALUES  ARE  FROM  FRAMES  SELECTED  VISUALLY  TO 
BE  PROPERLY  EXPOSED  BUT  DIFFICULT  TO  REDUCE. 

UNDERLINED  VALUES  OBTAINED  BY  SPECIAL  HAND  PROCESSING. 
OTHER  VALUES  OBTAINED  BY  SEMI -HANDS-OFF  METHOD. 

VALUES  IN  PARENTHESES  (1,9),  (15,6)  etc  REPRESENT  NUMBER  OF 
SEQUENTIAL  VALUES  OF  N (^)  AVERAGED  IN  A GROUP  FOLLOWED  BY 
NUMBER  OF  GROUPS  INVOLVED  IN  A PLOT. 
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3.  7 Each  value  of  a quoted  in  Table  2 represents  the  slope  of 

an  associated  plot  of  values  of  f (x  - xq)  vs  R (x-  xq).  All  such  plots 
are  presented  in  Appendix  A.  The  processed  data  from  each  horizon 
photograph  analyzed  (with  the  exception  of  the  one  failure)  is  presented 
in  the  appendix  as  three  plots  drawn  on  one  page.  These  are  the  plots  of 


N (<t>) 
N ( 


- N (*)  i _ 


f (x  - x ) vs  R (x  - x ). 
o o 


The  plots  are  for  the  blue,  the  green  and  the  red  image  of  the  horizon. 

The  plots  are  identified  by  date  and  Greenwich  mean  time.  The  letter  D in 
the  upper  right  hand  corner  identifies  a demonstration  run  and  the  letter  C 
a challenge.  The  values  of  total  scattering  coefficient,  which  is  the  slope 
of  the  plots  shown,  are  given  for  the  three  colors  as  are  the  associated 
aerosol  scattering  coefficients.  Certain  meteorological  information 
with  the  plots  and  other  supporting  information  is  given  in  Table  3. 


TABLE  3.  PHYSICAL  AND  METEOROLOGICAL  CONDITIONS 
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4.  DISCUSSION 


4.  1 The  demonstration  pictures  for  19  of  the  20  occasions  when 

/' 

the  horizon  was  photographed  give  reasonable  values  for  the  total  scat- 
tering coefficient  in  the  blue  the  green  and  the  red.  These  values  have 

t.  been  obtained  by  an  almost  automatic  process  of  data  reduction.  Operator 

involvement  occurs  in  the  selection  of  the  best  guess  position  of  the  horizon 
and  the  changes  of  this  guessed  position  in  the  iteration  process.  The 
twentieth  occasion.  May  18  at  0945,  required  hand  reduction  of  the  data 
and  yielded  unreasonable  results.  It  should  be  noted  that  the  challenge 
picture  selected  for  this  occasion  yielded  rather  reasonable  values  of 
the  scattering  coefficients. 


4-  2 The  photographs  selected  as  challenges  seem  to  produce 

reasonable  values  of  the  scattering  coefficients;  however,  there  is  only 
rather  general  agreement  between  simultaneous  challenge  and  the  demon- 
stration values.  The  values  on  May  24  at  1355  are  different  by  more  than 
an  order  of  magnitude.  In  this  case  the  demonstration,  data  reduction  was 
carried  out  the  automatic  process  and  involved  the  use  of  nine  adjacent 
pixels  for  data  processing.  In  this  case  the  values  of  scattering  coefficient 
have  been  determined  by  optical  effects  occuring  within  the  range  between 
approximately  1 and  4.  6 km  from  the  ship.  The  corresponding  challenge 
data  were  also  reduced  "automatically"  however,  using  the  first  nine 
pixels  below  the  horizon  did  not  work  well  in  the  "automatic"  process  so 
nine  sequential  groups  of  5 pixels  per  group  were  used.  In  this  case  the 
values  of  scattering  coefficient  have  been  determined  by  optical  effects 
occuring  within  the  approximate  range  0.  25  to  2.  4 km.  Inspection  of  the 


picture  used  for  the  challenge  shows  a dark  streak  on  the  sea  extending 
from  about  0.  16  km  from  the  ship  to  roughly  1.  3 km.  This  may  account 


for  the  results  obtained  with  this  film. 

4.  3 As  noted  before,  the  challenge  picture  for  31  May  at  1050 


could  not  be  analyzed  either  "automatically"  or  by  hand 


5.  CONCLUSIONS 

The  following  conclusions  seem  justified: 

1.  Color  photographs  of  the  sea  horizon  can  be  used  to  determine  scattering 
coefficients  at  X = 0. 44p,  X = 0.  55p,  X = 0.  64p. 

2.  The  semi-hands-off  calculations  demonstrate  that  completely  auto- 
matic data  analysis  is  feasible. 

3.  Such  an  analysis  system  could  be  used  in  a real  time  system  involving, 
for  instance,  scanning  or  a silicon  diode  array. 

4.  The  existing  films  could  be  used  as  a data  bank  for  simulations  of  the 
behavior  of  proposed  real  time  systems. 
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APPENDIX  A 


I 


Processed  Data 

In  the  following  pages  three  graphs  are  presented  for  each 
film  analyzed.  Each  page  is  specific  to  a particular  film  identified  by 
date  and  time.  In  addition  the  letter  D or  C identifies  demonstration 
pictures  or  challenge  pictures. 

The  graphs  are  of  the  function 

In  (N<0*>  - N<») 

( N(**) 

vs  R(£).  It  has  been  shown  that  if  the  graph  is  a straight  line  the  locus 
of  the  horizon  = o)  has  been  selected  correctly  and  the  slope  of  the 
straight  line  is  equal  to  the  scattering  coefficient. 

Twenty  pages  of  graphs  representing  the  reduced  data  for 
the  demonstration  films  are  presented  in  sequence  of  times  the  pictures 
were  taken. 

Following  this  nine  pages  of  challenge  graph  are  given.  Inter- 
spersed are  five  black  and  white  enlargements  of  the  color  pictures  used. 
The  demonstration  picture  included  is  for  3 June  at  1230.  Challenge 
pictures  are  18  May  at  0945,  21  May  at  1545,  22  May  at  1230  and  22  May 
at  1830.  Finally  a print  of  the  densitometer  trace  across  the  horizon  is 
given  for  the  blue  image  on  the  challenge  film  for  21  May  at  1545. 
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